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摘 要 : 风沙 运动 是 区 域 与 全 球 变化 研究 中 不 可 忽视 的 地 表 过 程 。 起 沙 阔 值 (临界 起 沙 风速 ,临界 起 沙 摩擦 速度 ) 
作为 判定 风沙 运动 是 否 发 生 的 关键 参数 ,是 风沙 运动 研究 的 核心 问题 。 利 用 塔克拉玛干 沙漠 腹地 野外 同步 观测 的 
风沙 运动 和 气象 资料 ,结合 数理 模型 计算 ,定量 评估 了 Stout, Kurosaki and Mikami (KM) 、 李 晓 风 和 张 宏 升 (LZ)、Mar- 
ticorena and Bergametti (MB ) „Shao and Lu( SL)5 种 起 沙 阔 值 判 定 方案 的 适用 性 ,并 利用 最 优 方案 确定 了 研究 区 新 的 
起 沙 闵 值 。 结 果 表 明 :(1) 5 种 方案 均 存 在 不 确定 性 ,KM 方案 确定 的 起 沙 阔 值 一 定 程 度 上 高 估 了 起 沙 时 长 与 沙 尘 
水 平 通 量 ,其 他 4 种 方案 则 相反 ,但 KM 方案 最 优 ,(2) 塔克拉玛干 沙漠 腹地 2 m 高 度 临 界 起 沙 风速 变化 范围 为 4.0~ 
6.0 ms ,临界 起 沙 摩擦 速度 变化 范围 为 0.24~0.36m:.s ;起 沙 阔 值 具 有 季节 性 差异 ,符合 夏季 > 秋季 > 春季 > 冬季 的 


变化 规律 。 
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风沙 运动 是 全 球 干 旱 - 半 干旱 区 重要 的 地 表 过 
程 之 一 。 风 沙 运动 可 以 引起 沙漠 化 与 土壤 风蚀 一， 
造成 道路 农田 与 村 庄 等 掩埋 **。 同 时 ,由 其 引起 
的 沙 尘 气 溶胶 排放 .传输 与 沉降 还 会 对 区 域 与 全 球 
气候 变化 生态 环境 "及 人 类 健康 "等 产生 
重大 影响 。 因 此 ,风沙 运动 的 相关 人 研究 受到 地 貌 . 气 
象 与 气候 .生态 与 环境 等 学 科 领 域 的 广泛 关注 P 

风 是 风沙 运动 发 生 的 重要 动力 因素 。 只 有 当 
风速 (w) 或 者 摩擦 速度 (w) 达 到 某 一 临界 值 时 ,风沙 
运动 才 会 发 生 , 该 值 即 临界 起 沙 风速 (w) 或 临界 起 
沙 摩 探 速度 (uv)。 临 界 起 沙 风 速 与 临界 起 沙 摩擦 速 
度 的 表征 意义 相同 ,前 者 可 通过 风速 和 风沙 运动 的 
观测 直接 确定 ,后 者 需 将 风速 换算 成 摩擦 速度 , 然 
后 配合 风沙 运动 观测 资料 确定 ;不同 学科 专业 的 
学 者 在 使 用 上 存在 差异 ,地 理学 科 的 学 者 多 使 用 临 
界 起 沙 风 速 ,而 气象 .环境 学 科 的 学 者 倾向 于 使 用 
临界 起 沙 摩擦 速度 ,本 研究 中 将 二 者 统称 为 起 沙 国 
值 。 起 沙 姜 值 不 仅 是 界定 风沙 运动 能 和 否 发 生 的 关 
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键 指标 ,也 是 土壤 风蚀 模型 和 沙尘暴 预报 模式 中 
的 核心 参数 ,决定 着 风蚀 强度 和 沙 尘 通 量 的 计算 
精度 …””"。 有 研究 表明 ,同等 条 件 下 起 沙 阔 值 降低 
32.1% 和 49.5% ,起 沙 量 可 增加 65.8% 和 69.296 * 。 
起 沙 阔 值 的 大 小 与 土壤 (如 土壤 湿度 .土壤 粒 径 和 
土壤 组 分 等 ) .植被 (如 植被 盖 度 等 ) 和 大 气 条 件 ( 如 
气温 和 空气 湿度 等 ) 等 密切 相关 …“ 后。 当 土 壤 湿 度 
大 于 对 起 沙 产生 影响 作用 的 临界 值 0.005 m -m° HT, 
起 沙 阅 值 随 土壤 湿度 增加 而 增 大 "1; 在 起 沙 参 数 化 
方案 中 ,起 沙 闵 值 被 表述 为 土壤 粒 径 的 函数 , 当 土 
壤 粒 径 为 74 um 左右 时 对 应 的 起 沙 阔 值 最 小 ,小 于 
或 者 大 于 该 粒 径 值 时 ,起 沙 阔 值 也 越 大 “ ;土壤 
中 因 含 有 盐分 而 形成 盐 壳 或 结 皮 可 使 起 沙 阔 值 增 
大 5; 植被 覆盖 可 减少 沙 源 供 给 ,降低 作用 于 地 表 
的 风 剪 切 力 , 进 而 降低 起 沙发 生 几 率 , 增 大 起 沙 阔 
值 ““* ;空气 湿度 的 增 大 可 提高 表层 土壤 湿度 或 者 
增加 土壤 颗粒 间 的 粘 清 度 ,进而 导致 起 沙 阔 值 升 
pj: 气温 的 变化 可 以 改变 空气 湿度 ,进而 影响 起 
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沙 阔 值 汪 3 。 起 沙 阔 值 的 判定 方法 已 开展 了 较 多 研 
究 , 李 晓 岚 站 将 其 分 为 观测 法 包括 风 洞 实验 "2 沁 和 
外 场 试验 ”统计 法 -及 参数 化 "I。 各 种 方 
法 对 比 研 究 的 相关 工作 开展 尚 少 ,因此 , 哪 种 方法 
更 有 优势 至 今 尚 不 明确 。 

我 国有 荒漠 化 土地 261.16x10+ km? ,是 受 风 沙 运 
动 危害 较为 严重 的 国家 之 一 ?5。 塔 克拉 玛 干 沙漠 
为 我 国 第 一 大 沙漠 ,风沙 活动 强度 在 世界 上 首届 一 
指 汪 。 同 时 ,其 所 处 的 塔里木 盆地 是 新 疆 社 会 .经 
济 文化 等 的 重要 承载 区 ,塔克拉玛干 沙漠 及 周边 
地 区 分 布 有 新 疆 53.5% 的 绿洲 ,居住 着 新 疆 48.07% 
WAO , 储 茂 有 丰富 的 油气 等 矿产 资源 小。 因此 ， 
开展 该 区 域 起 沙 阔 值 研究 ,为 风沙 灾害 及 沙 尘 天 气 
发 生 提供 预报 预警 ,对 当地 的 经 济 社会 发 展 和 生态 
环境 保护 具有 重要 的 意义 。 

关于 塔克拉玛干 沙漠 起 沙 阔 值 方面 的 研究 工 
WE ,主要 表现 在 如 下 几 个 方面 : 陈 渭南 等 在 沙漠 
北 缘 肖 塘 地 区 通过 人 工 目 测 风沙 运动 发 生 情况 , 配 
合同 步 观测 的 风速 数据 ,确定 2 m 高 度 瞬 时 冲击 临 
界 起 沙 风 速 为 5.0 ms"', 该 值 被 广泛 应 用 于 塔 克 拉 
玛 干 沙漠 风沙 运动 研究 ;Ishizuka 等 * 在 沙漠 南 缘 
策 勒 地 区 通过 光学 传感器 测 得 3.8 m 高 度 临 界 起 沙 
风速 分 别 为 7.5 mes CEP) ffll9.5 m-s Gh) ; Ku- 
rosaki 等 ”基于 沙漠 周边 气象 站 风速 和 沙尘暴 资 
料 ,统计 得 到 该 区 域 10m 高 度 临 界 起 沙 风速 为 5.2~ 
8.2 ms ;Yang 等 ”利用 沙漠 腹地 所 测 沙 粒 跃 移 数 
据 , 计 算得 到 2 m 高 度 起 沙 风 速 为 3.5~10.9 ms; 
Zhou 等 "基于 外 场 观测 资料 对 比 了 几 种 常用 起 沙 
闵 值 判定 方法 的 差异 。 尽 管 如 此 ,塔克拉玛干 沙漠 
起 沙 阔 值 的 研究 仍 存在 以 下 不 足 之 处 :观测 期 较 
短 ,缺少 对 起 沙 阔 值 季节 变化 的 认识 ;起 沙 阔 值 判 
定 方 法 在 研究 区 的 适用 性 未 进行 评估 。 本 研究 基 
于 塔克拉玛干 沙漠 腹地 野外 长 期 观测 数据 ,评估 5 
种 常用 起 沙 阔 值 判定 方法 的 适用 性 , 遵 选 出 最 优 判 
定 方 法 ,进而 确定 新 的 起 沙 阔 值 ,为 研究 区 风沙 运 
动 的 判定 提供 更 精准 的 判 据 。 


1 研究 区 概况 


观测 试验 在 塔克拉玛干 沙漠 腹地 塔 中 地 区 开 
JÈ (39°00' N, 83°38' E, 15:1 1103 m; 图 1)。 塔 中 地 
区 下 垫 面 均 为 流沙 地 表 ; 地 表土 壤 粒 径 集 中 分 布 在 


63-250 hm 之 间 , 中 值 粒 径 为 147 um“. XE EC 
34 55 25 2 r8] vb 28 55 28 [8] HR [8] 7p 48 , vb 2E 3E EON 
NNE-SSW zX EN-SW 3E [i] ,沙化 与 垄 间 地 的 相对 高 
度 约 40~50 m, Æ [R] Jb 952^] 1-3 km, 长 约 2~5 km, H 
ARX PE, KARRE er T 2 2& BECK 
沙 垄 间 的 平坦 沙 地 上 ,利用 塔 中 气象 站 1996—2017 
年 气象 观测 资料 ,得 到 研究 区 年 均 气温 为 11.7 C, 
有 观测 记录 以 来 最 高 气温 为 46.0 ,最 低 气温 为 
-32.6% ;年 均 降 水 量 为 27.6 mm,3 一 8 月 的 降水 约 
di EERE HI 87.696; 4.35] 28 A EIK 3741.8 mm; 年 
15] GS JJ 2.2 ms ,年 均 大 风 日 数 为 10.7 d, FER 
中 在 春 、 夏 季 ; 年 均 沙尘暴 . 扬 沙 和 浮尘 日 数 分 别 为 
17.0 d,68.0 4 和 122.0 d, 春 、 夏 季 为 高 发 期 。 


2 数据 与 方法 


2.1 风沙 运动 数据 

借助 Sensit 压 电 式 风 刨 传感器 和 Big Spring 
Number Eight (BSNE ) 集 沙 仪 开 展 风沙 运动 观测 (图 
1)。 当 沙 粒 撞击 风蚀 传感器 的 晶片 时 ,会 输出 脉冲 
信号 给 数据 采集 器 (Campbell CR1000) ,记录 下 撞击 
颗粒 的 数量 ;无 沙 粒 撞击 时 则 记录 为 0。 数 据 采 集 
频率 设 为 1 Hz, 可 提供 1s 时 间 步 长 撞击 颗粒 数 的 连 
续 观 测 数 据 。 当 某 1 s 的 撞击 颗粒 数 大 于 0 时 , 则 记 
录 为 起 沙 时 长 1 s, 进 而 可 统计 每 次 沙尘暴 天 气 过 程 
的 总 撞击 颗粒 数 和 起 沙 时 长 。 风 蚀 传感器 安装 高 
度 为 5 cm 和 10 cm, 本 研究 仅 使 用 5 em 高 度 的 测量 
数据 。Sensit 风蚀 传感器 经 过 长 期 野外 试验 验证 ， 
对 人 研究 区 起 沙 过 程 的 测量 效果 良好 四。 为 了 弥补 
风蚀 传感器 不 能 测量 沙 尘 输送 通 量 的 缺陷 ,在 其 一 
侧 200~300 cm 处 安装 BSNE 梯度 集 沙 塔 , 且 2 组 仪 
器 的 连 线 与 研究 区 主 风 向 垂直 。BSNE 集 沙 塔 的 高 
度 为 200 cm ,安装 有 6 组 集 沙 仪 , 集 沙 盒 进 沙 口中 间 
距离 地 面 分 别 为 S cm 10 cm .20 cm 、50 cm 、100 cm 
fll 200 cm; 采 样 频率 为 单 次 沙尘暴 或 扬 沙 天 气 过 
程 。2 组 仪器 测量 数据 相互 配合 可 以 获取 观测 点 
100 cmx200 cm 截面 上 不 同时 间 步 长 的 沙 侍 水 平 通 
量 , 用 以 验证 沙 尘 水 平 通 量 参数 化 方案 的 计算 值 ”。 
2.2 近 地 层 微 气象 和 沙 侍 浓度 数据 

近 地 层 微 气象 借助 10 m 梯 度 探测 平台 开展 观测 
(图 1)。 该 平台 架设 有 0.5 m1m.2m.4m 和 10m5 
层 梯度 风速 (Metone 010C) 、 风 向 (Metone 020C) 和 


注 :(a) 试 验 点 塔 中 位 置 ;(b) 试 验 场 全 景 ;(c)Sensit 压 


电 式 风蚀 传感器 ; (d) Big Spring Number Eight(BSNE ) 集 沙 仪 。 


图 1 塔克拉玛干 沙漠 试验 点 位 置 及 观测 仪器 


Fig. 1 Location of experimental site in the Taklimakan Desert and the instruments that were used in the field experiments 


空气 温 湿 度 (Vaisala HMP45C ) 传 感 器 ,并 在 3 m Sj 
度 架设 了 沙 尘 浓度 测量 仪 (Grimml.108)。 同 时 该 
平台 还 设 有 1.5 em、5 cm 和 10 cm 3 层 土壤 湿度 探头 
(Campbell CS616) 。 微 气象 测量 数据 统一 存储 在 
CR1000 型 数据 采集 需 中 (Campbell) ,可 提供 1s、 
1 min、 30 min 和 1p 时 间 步 长 的 测量 数据 。 借 助 
Grimm1.108 开展 沙 尘 浓度 的 测量 ,该 仪器 可 测量 
20 um 以 下 沙 尘 粒 子 的 浓度 ,测量 频率 为 6s, 然 后 平 
1573 1 min fll 5 min 的 值 。 沙 侍 天 气 数据 源 于 塔 中 气 
象 站 2008 年 2 月 一 2018 年 3 月 天 气 现象 观测 记录 。 
2.3 起 沙 阐 值 判定 方法 

基于 试验 内 容 , 选 取 Stout2 方案 代表 野外 观测 
法 ;Kurosaki and Mikami (KM) FE , ZEE Jt FUE 
JF (OZ) 7 Jr RAR BE YE IE ; Marticorena and Ber- 
gametti (MB) 方案 、Shao and Lu(SL) * 7; 364 
数 化 。 具 体 如 下 : 

CL) 观测 法 

Stout ^ 基于 Sensit 风蚀 传感器 测量 的 风沙 运动 


数据 ,发 展 了 起 沙 风 速 的 判定 方法 : 
u, 7 üg d (y) (1) 
式 中 :为 分 钟 临 界 起 沙 风速 ms ) ;五 为 分 钟 平均 
风速 (ms"');o 为 分 钟 平均 风速 的 标准 差 ;y 为 风沙 
运动 强度 系数 ,用 单位 时 间 内 起 沙发 生 时 间 所 占 比 
表示 , 取 值 范围 在 0~1 之 间 , 当 y=1 时 ,表示 风沙 运 
动 持 续 发 生 ;y=0 时 ,表示 无 风沙 运动 发 生 ; PNY 
的 正 态 分 布 孔 数 。 该 方法 已 被 用 于 多 个 地 区 起 沙 
阔 值 的 判定 2 。 
(2) 统计 法 
Kurosaki and Mikami ^" 基于 风速 和 沙尘暴 观测 
数据 ,发 展 了 起 沙 阔 值 的 判定 方法 : 
P,=n,/N, X100 (2) 
式 中 :P; 为 百分比 ;nj 为 风速 等 级 为 i 时 风沙 运动 发 
生 的 频次 ;NN 为 风速 等 级 为 i 的 总 频次 ;风速 等 级 i 的 
分 组 间隔 为 0.2m…s';P; 为 50% 时 的 风速 等 级 定义 
为 临界 起 沙 风 速 w。 
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李 晓 风 和 张 宏 升 基于 内 蒙古 科尔沁 沙 地 起 
沙 试验 数据 ,将 一 定 观 测 高 度 沙 侍 浓度 开始 迅速 且 
连续 增 大 时 对 应 的 风速 判定 为 临界 起 沙 风 速 。 

(3) 参数 化 

Marticorena and Bergametti * , Shao and Lu ^ 发 
展 的 起 沙 阔 值 参 数 化 方案 是 当前 沙尘暴 模式 中 使 
用 最 为 广泛 的 方案 。 方 案 均 可 用 下 式 表达 : 


wu. (d, À, w) = us (d) f (A) f. (w) (3) 
在 MB 方案 中 ， 
0.129K 0.08 « B «10 


us (d)- [0.06178 - 10] 
0.129k (1 - 0.08586" ) , B>10 
(4) 
PORE = (5) 
In[0.702255cnz,)"*] 
1 , w«w' 
fw) ES ln 十 1.21(w E w! Maig ! w > w' (6) 


式 中 :ws(g) 为 光滑 地 表 临 界 起 沙 摩擦 速度 ;A(A) 为 地 
表 粗 糙 元 修正 方程 ;Aiko) 为 土壤 湿度 修正 方程 。B= 
ad' * b, aa 和 x* 分 别 为 1331 cm“*、0.38 和 1.56; K= 
[(cgd/p)“(1+0.006/og(d)”)]”,g 为 重力 加 速度 , 取 值 
9.81 ms ”7,0 为 沙 粒 密度 , 取 值 2650 kg*m i,p 为 空气 
密度 , 取 值 1.02 kgmi,d 为 沙 粒 粒 径 , 取 值 147 um; 
为 为 地 表 粗 糙 度 , 取 值 0.2478 em; zo WIGI He RA 
度 , 取 值 0.01372 cm;w 为 土壤 重量 含水 量 ;w' 为 土壤 
重量 含水 量 是 否 影 响起 沙 的 临界 值 [w'=0.0014(% 
clay)^0.17(96clay) | * 


在 SL 方案 中 ， 
oed I0.5 

nores] D 

f= -mo AU +m, B, Ay? (8) 


e?" : w « 0.03 

/= pm f w > 0.03 

XP : Ave mo 和 B. 为 系数 , 取 值 分 别 为 0.0123、 

0.165 g*s7,0.5,1 1190; A--0.35In(1-a) , a 为 植被 覆 
盖 指 数 , 取 值 0;w 为 土壤 体积 含水 率 。 

由 于 2 组 方案 计算 的 均 为 临界 起 沙 摩擦 速度 

ua, 可 利用 下 式 将 其 换算 为 ,以 方便 与 其 他 方法 确 

定 的 起 沙 阔 值 进行 比较 : 


ME 
LET (2) (10) 


CB Uc EAR ,BUBLO.4;z yu s BE UE 2.0 mo 
2.4 沙 尘 水 平 通 量 确定 方法 
沙 尘 水 平 通 量 利用 Owen 方案 计算 : 


Q;(d)= Z l 3i (11) 


式 中 : 0,(q) 为 粒 径 为 ;的 粒 径 组 沙 尘 水 平 通 量 (kg， 
m'*min ^) ,与 粒 径 组 ;的 起 沙 阔 值 对 应 ,由 于 笔者 试 
验 获 取 的 起 沙 阔 值 代 表 了 整体 粒 径 组 ,因此 QO(dg) 即 
为 0;E 为 可 发 生起 沙 的 地 表 比 重 ,研究 区 取 值 1;c 
为 系数 , 取 值 0.8。 


3 结果 与 分 析 


3.4 起 沙 阅 值 判 定 方法 适用 性 评估 

根据 试验 数据 的 完整 性 ,本 人 研究 选择 2009 4E. 7 
月 17 日 的 沙尘暴 过 程 为 例 ,验证 各 起 沙 阔 值 判定 方 
法 在 研究 区 的 适用 性 ,图 2 给 出 了 本 次 沙尘暴 过 程 
风速 与 风沙 运动 变化 情况 ,由 图 2 可 知 , 本 次 沙尘暴 
过 程 中 风速 从 4:00 一 5:00 左 右 开 始 增 大 ,在 6:00 左 
右 达 到 了 5.0 ms ,此 时 已 可 监测 到 较 弱 的 起 沙 活 
动 ;风速 在 9:00 左 右 增 加 到 了 6.0 ms 以 上 ,风沙 运 
动 强度 迅速 增 大 ,持续 至 19:00 左 右 , 随 着 风速 的 降 
AI ,风沙 运动 减弱 至 结束 。 本 次 沙尘暴 过 程 共 发 生 
风沙 运动 735 min, 对 应 的 最 小 分 钟 平 均 风 速 为 2.5 
ms ,最 大 为 11.2 ms ,均值 为 7.3 meso 

图 3 分 别 给 出 了 Stout KM LZ .MB 和 SL 方 案 确 
定 的 临界 起 沙 风 速 ,Stout 方 案 确 定 的 分 钟 临界 起 沙 
风速 的 变化 范围 为 3.1~9.1 ms ,均值 为 6.4 m:s 5 
分 钟 临界 起 沙 摩 擦 速 度 的 变化 范围 为 0.19~0.54 m+ 
s ,均值 为 0.38 ms'; 以 分 钟 临界 起 沙 风速 为 风沙 
运动 判 据 ,本 次 沙尘暴 过 程 风 沙 运 动 发 生 时 长 203 


60 - 12 
um 起 沙 时 间 

E 一 风速 m 
E 40 $& «m 
[2i T 
ES E 
8 R 
La 20 4 R 

0 0 

00:00 04:00 08:00 12:00 16:00 20:00 00:00 


时 刻 


图 2 2009 年 7 月 17 日 沙尘暴 过 程 风 沙 运动 与 风速 变化 


Fig.2 Changes of wind-blown sand movement and wind 


speeds during sandstorm process on July 17, 2009 
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图 3 2009 年 7 月 17 日 沙尘暴 过 程 Stout .KM LZ. MB 和 SL 方案 确定 的 临界 起 沙 风 速 
Fig. 3 Threshold velocities for Stout, KM, LZ, MB and SL schemes during sandstorm process on July 17, 2009 


min, 为 实际 发 生 时 长 的 27.6% ;以 均值 为 判 握 ,风沙 
运动 发 生 时 长 517 min, 为 实际 发 生 时 长 的 70.3%。 
根据 KM 的 判定 方法 ,本 次 沙尘暴 过 程 中 临界 起 沙 
风速 约 为 4.8 ms ,临界 起 沙 摩擦 速度 约 为 0.29 me 
s ;以 该 值 为 风沙 运动 判 据 ,本 次 沙尘暴 过 程 风沙 运 
动 发 生 时 长 756 min ,为 实际 发 生 时 长 的 102.9% 。 
LZ 方 案 判定 的 临界 起 沙 风 速 约 为 6.8 ms ,临界 起 
沙 摩 氛 速 度 约 为 0.41 ms'!; 以 该 值 为 判 据 , 本 次 沙 
侍 暴 过 程 风沙 运动 发 生 时 长 457 min ,为 实际 发 生 时 
长 的 62.2%。 由 于 研究 区 下 热 面 和 土壤 湿度 变化 极 
小 ,MB 方案 计算 的 分 钟 临界 起 沙 风 速 几乎 无 变化 ， 
值 约 为 6.6 ms ,临界 起 沙 摩 擦 速度 约 为 0.39 mes; 
以 该 值 为 风沙 运动 判 据 , 本 次 沙尘暴 过 程 风沙 运动 
发 生 时 长 495 min ,为 实际 发 生 时 长 的 67.3%。SL 方 
案 计 算 的 分 钟 临界 起 沙 风 速 变化 范围 为 5.4~6.0 mm: 
s ,均值 为 5.6 ms ;分 钟 临界 起 沙 摩擦 速度 变化 范 
围 为 0.32~0.36 ms ,均值 为 0.33 m: s ;以 分 钟 临界 
起 沙 风 速 为 判 据 , 本 次 沙尘暴 过 程 风 沙 运动 发 生 时 


长 614 min; 以 均值 为 判 据 , 风 沙 运动 发 生 时 长 594 
min, 分 别 为 实际 发 生起 沙 时 长 的 83.5% 和 80.8% 。 
上 述 结果 表明 ,除了 KM 方案 ,其 余 4 种 方案 均 不 同 
程度 地 低估 了 风沙 运动 时 长 。 

为 了 进一步 验证 各 起 沙 阔 值 判定 方案 在 研究 
区 的 适用 性 ,本 研究 基于 各 方案 确定 的 临界 起 沙 阔 
1E (Stout &II SL 方案 为 均值 ) ,结合 Owen55 沙 侍 水 平 
通 量 计算 方案 ,计算 了 2009 年 7 月 17 日 沙尘暴 过 程 
分 钟 沙 尘 水 平 通 量 ,并 与 实测 结果 进行 了 对 比 ( 图 
4)。 本 次 沙尘暴 过 程 实测 沙 尘 水 平 通 量 为 271.0 
kg.m- ,基于 Stout KM , LZ , MB 和 SL 方案 确定 的 起 
沙 阔 值 计算 的 沙 尘 水 平 通 量 分 别 为 178.3 kgs m^, 
314.1 kg*m'',122.9 kg*m'',179.2 kg*m 和 244.7 kg- 
nm 分 别 为 实测 值 的 65.8% , 115.896 ,45.396 , 66.196 
和 90.2%。 与 风沙 运动 时 长 相似 ,基于 KM 方案 确定 
的 起 沙 阔 值 计算 的 沙 尘 水 平 通 量 在 一 定 程度 上 高 
估 了 本 次 沙尘暴 过 程 中 风沙 运动 强度 ,其余 4 种 方 
案 均 不 同 程度 地 低估 了 本 次 沙尘暴 过 程 中 风沙 运 
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图 4 2009 年 7 月 17 日 沙尘暴 天 气 过 程 观测 与 计算 的 沙 尘 水 平 通 量 比较 


Fig. 4 Comparisons of observed and calculated horizontal dust fluxes with Stout, KM, LZ, MB and SL schemes during sandstorm 


process on July 17, 2009 


动 强度 。 相 关 性 分 析 表 明 ,KM 方案 效果 最 优 ,相关 
系数 7 为 0.90, 依 次 为 9L 方 案 (r=0.87) Stout J Z (r= 
0.86) .MB 方案 (r=0.82) .LZ 方案 (r=0.80)。 研 究 表 
BH ,脉动 风速 是 导致 沙 粒 运动 发 生 的 重要 动力 机 
制 ,上 且 脉动 风速 与 起 沙 量具 有 较 好 的 一 致 性 “。 本 
研究 中 使 用 的 风速 数据 时 间 步 长 为 1 min, 这 在 一 定 
程度 上 会 屏蔽 掉 部 分 风速 脉动 信息 ,进而 导致 估算 
的 起 沙 时 长 与 沙 尘 水 平 通 量 与 实际 情况 存在 一 定 
差异 。 利 用 上 述 5 种 起 沙 阔 值 判定 方案 ,进一步 统 
计 计 算 了 2009 年 研究 区 其 他 8 次 沙尘暴 天 气 过 程 
起 沙 时 间 与 沙 侍 水 平 通 量 ( 表 1 和 表 2) ,结果 与 
2009 年 7 月 17 日 沙尘暴 过 程 一 致 。 综 上 所 述 ,KM 
方案 最 为 适合 研究 区 起 沙 阔 值 的 判定 。 由 于 土壤 
和 植被 等 下 垫 面条 件 存 在 区 域 差异 ,KM 方案 在 其 
他 下 垫 面 (如 荒漠 .戈壁 . 干 酒 河 / 湖 床 等 ) 是 否 仍 然 


为 起 沙 阔 值 的 最 优 判定 方案 有 待 验 证 ,因此 ,未 来 
有 必要 开展 不 同 下 垫 面 .不同 起 沙 闵 值 判定 方案 的 
对 比 研究 工作 。 
32 塔克拉玛干 沙漠 起 沙 阅 值 

利用 KM 方法 确定 了 人 研究 区 2008 年 3 月 一 2018 
年 2 月 101 次 沙 侍 天 气 过 程 的 起 沙 病 值 (图 5), 结 
表明 ;2 m 高 度 临界 起 沙 风 速 变 化 范围 为 4.0~6.0 m- 
s ,均值 为 5.0 ms ;临界 起 沙 摩 擦 速度 变化 范围 为 
0.24~0.36 m" s” ,均值 为 0.30 m- s^; 临界 起 沙 风速 分 
布 较为 集中 ,73.2% 的 值 位 于 4.6~5.2 ms zz] oe 
沙 阔 值 具有 显著 的 季节 变化 , 春 . 夏 .秋冬 季 的 临 
界 起 沙 风 速 分 别 为 4.9 m-s',5.4 m:s',5.1 m:s 和 
4.6 m*s ,临界 起 沙 摩擦 速度 分 别 为 0.29 ms ! 0.32 
m:s/,030 m:s ! 410.27 m.s ,符合 夏季 > 秋季 > 春 
季 > 冬 季 的 变化 规律 ( 表 3)。 前 人 的 研究 表明 ,土壤 
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表 1 塔克拉玛干 沙漠 不 同 沙尘暴 天 气 过 程 观测 与 计算 的 起 沙 时 间 比 较 


Tab.1 Occurrence times of sand movement from observed and calculated with emission thresholds determining methods 


during different sandstorms in the Taklimakan Desert 


时 间 / 年 -月 -日 观测 值 /min Stout 方 案 /min KM 方案 /min ILZ 方 案 /min MB 方案 /min SL 方案 /min 
2009-03-10 788 680 833 631 626 719 
2009-03-26 663 205 731 91 29 238 
2009-04-16 501 404 523 368 363 428 
2009-04-29 872 756 907 678 786 806 
2009-05-20 914 633 965 532 685 860 
2009-06-05 589 238 619 152 238 507 
2009-06-30 870 707 885 699 743 850 
2009-08-05 228 115 300 95 100 168 


表 2 塔克拉玛干 沙漠 不 同 沙 尘 暴 天 气 过 程 观测 与 计算 的 沙 尘 水 平 通 量 比较 


Tab. 2 Horizontal dust fluxes from observed and calculated with emission thresholds determining methods during different 


sandstorms in the Taklimakan Desert 


时 间 / 年 -月 -日 ”观测 值 /(kg'm) ”Stout 方案 /(kg:m') KM 方案 /(kg*m')” ”12 方案 /(kg*m') ”MB 方案 /(kg:m') ”SL 方案 /(kg*m') 
2009-03-10 416.4 311.2 442.8 273.4 263.9 357.3 
2009-03-26 45.1 114 60.3 5.5 23 134 
2009-04-16 285.6 232.6 3513/7 209.6 203.9 250.8 
2009-04-29 368.5 2844 379.3 237.4 308.2 324.6 
2009-05-20 162.8 107.5 183.9 90.9 116.4 144.3 
2009-06-05 81.9 25.9 100.1 11.4 25.9 68.7 
2009-06-30 412.9 328.8 426.3 316.5 365.3 382.3 
2009-08-05 46.2 21.7 49.3 16.2 17.5 32.8 

120 气 湿度 的 季节 差异 较 小 ,导致 该 区 域 起 沙 阔 值 季节 
100 变化 较 微 弱 。Kurosaki and Mikami” 利用 塔 克 拉 玛 

干 沙漠 周边 气象 站 1988—2005 年 风速 和 沙尘暴 观 

B 测 数据 ,统计 得 到 该 区 域 10 m 高 度 临 界 起 沙 风 速 的 
m 60 变化 范围 为 5.2~8.2 m-s”, HHEN 6.7 m* s^ ;本 研究 
ES m 的 临界 起 沙 风 速 换 算 成 10 m 高 度 为 5.0~7.6 m-s”, 
均值 为 6.3 m- s^, MAIR F Kurosaki and Mikami 的 统 

20 计 结 果 。 分 析 其 原因 ,Kurosaki and Mikami 使 用 的 

1 数据 为 沙尘暴 时 时 风速 ,而 沙尘暴 是 由 强烈 的 风沙 


39 41 43 45 4.7 49 51 53 5.5 57 59 61 
临界 起 沙 风速 /m.s2) 


图 5 研究 区 临界 起 沙 风 速 频率 累积 分 布 


Fig. 5 Cumulative distribution of threshold velocities in 


study area 


湿度 .空气 湿度 和 气温 的 升 高 均 能 使 起 沙 阔 值 增 
KOS, WK EREE ,水汽 压 和 气温 的 季节 变 
化 与 起 沙 阔 值 一 致 ,使 得 起 沙 阔 值 呈现 上 述 季 闻 变 
化 。 由 于 塔克拉玛干 沙漠 降水 稀少 ,土壤 混和 度 与 空 


运动 引起 的 ,其 对 应 的 风速 必然 高 于 临界 起 沙 风 
速 ;本 研究 的 结果 是 基于 5 em 高度 风 沙 运 动 数据 与 
同步 的 风速 数据 统计 得 出 ,可 更 加 真实 地 反映 研究 
区 起 沙 阔 值 的 变化 情况 。 


4 讨论 


本 研究 基于 野外 观测 数据 评 佑 了 5 种 常用 起 沙 
闵 值 判定 方法 ,虽然 上 述 5 种 方法 均 在 世界 不 同 地 
区 的 风沙 运动 研究 中 得 以 应 用 ,然而 将 各 方法 用 于 
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表 3 塔克拉玛干 沙漠 起 沙 阅 值 季节 变化 


Tab.3 Seasonal variations of threshold velocities and friction velocities in the Taklimakan Desert 


季节 平均 临界 起 沙 风速 ”平均 临界 起 沙 摩擦 速度 “平均 土壤 湿度 平均 大 气 水 汽 压 平均 气温 
K(m*s" /(m*s) /(m'*m?) /kPa PC 

春季 4.9 0.29 0.018 33 174 

夏季 5.4 0.32 0.021 6.7 28.2 

秋季 54 0.30 0.019 5.6 18.7 

冬季 4.6 027 0.017 1.6 -1.8 

表 4 不 同 沙漠 / 沙 地 起 沙 阔 值 比较 
Tab. 4 Comparisons of threshold velocities and friction velocities of different deserts and sands 
m 观测 高 度 临界 起 沙 风速 临界 起 沙 摩擦 速度 年 均 降水 植被 覆盖 文献 来 源 
/m /ms ) Im» s~) /mm 1% 

塔克拉玛干 沙漠 2 4.0~6.0 0.24-0.36 25.7 0.0 本 研究 
塔克拉玛干 沙漠 10 5.0~7.6 0.30~0.45 25.7 0.0 本 研究 
古 尔 班 通 古 特 沙漠 2 11.1-13.9 - 156.0 32.3 47] 
科尔沁 沙 地 4 6.9~11.5 0.55~0.70 366.0 50.0 24] 
Mescalero Sands, USA 2 10.0-11.9 - 300.0 - 48] 
浑 善 达 克 沙 地 4 10.0 0.6 365.0 20.0 49] 
敦煌 沙 地 10 7.0 0.5 23.2 0.0 20] 
戈壁 沙漠 10 11.8~15.8 一 一 一 31] 


同一 区 域 起 沙 阔 值 的 判定 工作 开展 较 少 。 通 过 比 
较 ,发 现 各 方法 之 间 结 果 的 差异 性 较为 显著 ,为 风 
沙 运动 的 判定 带 来 较 大 不 确定 性 ,因此 开展 各 方法 
的 适用 性 评估 具有 必要 性 。 本 研究 中 ,Stout 方 案 利 
用 野外 风沙 运动 资料 可 以 快速 获取 起 沙 阔 值 ,但 该 
方法 是 基于 风速 脉动 的 正 态 分 布 建立 的 ,在 真实 时 
外 状态 下 风速 脉动 的 分 布 是 否 符合 正 态 分 布 有 待 
商 榨 ;KM 方案 选择 将 某 一 风速 发 生 风 沙 运 动 的 概 
率 达 到 50% 时 定 为 临界 起 沙 风速 ,从 统计 学 的 角度 
具有 一 定 的 合理 性 ,但 是 该 方案 需要 大 量 样 本 数 
据 ;LZ 方 案 基 于 沙 尘 浓度 变化 (时 间 ) 来 判定 起 沙 阔 
值 ,具有 一 定 的 主观 性 ;MB 和 LS 方案 的 理论 基础 较 
为 完善 ,但 计算 过 程 较 为 繁琐 ,需要 获取 的 参数 较 
多 ,日 很 多 系数 源 于 风 洞 实验 或 外 场 的 观测 试验 ， 
系数 和 参数 均 存 在 地 域 适用 性 问题 。 因 此 , 至今 仍 
缺乏 一 种 快速 .简便 ,准确 的 起 沙 羡 值 判定 方法 。 
与 国内 外 不 同 地 区 沙漠 的 起 沙 阔 值 比较 可 知 ， 
塔克拉玛干 沙漠 的 起 沙 阔 值 最 小 ( 表 4) 2 7” ,日 
变化 范围 较 罕 。 虽 然 降水 .植被 覆盖 条 件 与 敦 烛 沙 
地 类 似 , 但 仍 略 低 于 敦煌 沙 地 7.0 ms 和 0.5 m-s” 
的 起 沙 阔 值 。 与 其 他 地 区 的 沙漠 相 比 ,塔克拉玛干 
沙漠 降水 量 较 小 ,气候 干燥 , 旦 无 植被 覆盖 ,为 风沙 


运动 的 发 生 提 供 了 有 利 的 气候 与 下 扑面 条 件 ; 塔 元 
拉 玛 干 沙漠 地 表土 壤 粒 径 以 细 砂 和 极 细 砂 为 主 , 平 
均 粒 径 在 100 um 左右 , 正 处 于 起 沙 阔 值 的 低 值 
pc ,为 起 沙 提供 了 有 利 的 土壤 条 件 , 上 述 因 素 
是 塔克拉玛干 沙漠 起 沙 阔 值 较 小 的 形成 原因 。 


5 结论 

通过 野外 观测 获取 的 风沙 运动 时 长 和 沙 侍 水 
平 通 量 数据 评估 了 5 种 常用 起 沙 阔 值 判定 方案 在 塔 
克拉 玛 干 沙漠 腹地 的 适用 性 ,并 基于 最 优 判定 方案 
确定 了 塔克拉玛干 沙漠 腹地 的 起 沙 阔 值 。 具 体 结 
论 如 下 : 

(1) Stout, LZ .MB 和 SL4 组 起 沙 阔 值 判定 方案 
均 不 同 程度 高 估 了 研究 区 起 沙 阔 值 ,使 用 上 述 起 沙 
国 值 确定 研究 区 风沙 运动 时 长 和 沙 尘 水 平 通 量 时 
将 存在 不 同 程度 的 低估 现象 ;基于 KM 方案 判定 的 
起 沙 阔 值 确定 的 风沙 运动 时 长 与 实测 值 最 为 接近 ， 
沙 尘 水 平 通 量 与 实测 值 接近 且 相 关 性 最 好 。 由 此 
可 以 认定 KM 方案 为 判定 塔克拉玛干 沙漠 腹地 起 沙 
阔 值 的 最 优 方 案 。 

(2) 利用 KM 方案 统计 得 到 塔克拉玛干 沙漠 腹 
地 2m 高 度 临 界 起 沙 风速 变化 范围 为 4.0~6.0 m-s”, 
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均值 为 5.0 ms ;临界 起 沙 摩擦 速度 变化 范围 为 
0.24~0.36 ms ,均值 为 0.30 ms ;起 沙 阔 值 符合 
季 > 秋 季 > 春 季 > 冬 季 的 变化 规律 。 
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Comparative study of the applicability of emission threshold 
determining method in the Taklimakan Desert 


YANG Xinghua", MA Mingjie, ZHOU Chenglong, HE Qing’ 
(1. School of Geographical Sciences, Shanxi Normal University, Taiyuan 030032, Shanxi, China; 
2. National Observation and Research Station of Desert Meteorology, Taklimakan Desert in Xinjiang, 


Urumqi 830002, Xinjiang, China) 


Abstract: Aeolian sand movement is a crucial surface process in the study of regional and global changes. As a 
key parameter in determining wind-sand movement occurrence, the emission threshold (threshold velocity, threshold 
friction velocity) is a core issue in the study of wind-sand movement. Taking the Taklimakan Desert as the target 
area, using synchronously observed wind- sand movement and meteorological data in the field, combined with 
mathematical model calculations, the applicability of five emission threshold judgment schemes from Stout, 
Kurosaki and Mikami (KM), Li Xiaolan and Zhang Hongsheng (LZ), Marticorena and Bergametti (MB), and 
Shao and Lu (SL) was quantitatively assessed, and the optimal scheme was used to determine the new emission 
threshold of the Taklimakan Desert. The following results are presented. (1) Uncertainties exist in all five 
schemes. The emission threshold determined by the KM scheme overestimates the dust emission time and the 
dust level flux to a certain extent. The four other schemes are opposite, but the KM scheme is the best. (2) The 2 
m threshold velocity in the hinterland of the Taklimakan Desert varies from 4.0 m*s ' to 6.0 ms ', and the 
threshold friction velocity varies from 0.24 m-*s ' to 0.36 m-s '. Moreover, the emission threshold has seasonal 
differences, which is consistent with summer > autumn > spring > winter. 

Keywords: wind-blown sand movement; threshold velocity; threshold friction velocity; Taklimakan Desert 
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